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ABSTRACT 

We investigate the environmental dependence of the mass-metallicity relation at z = 2 with MOSFIRE/Keck 
as part of the ZFIRE survey. Here, we present the chemical abundance of a Virgo-like progenitor at z = 2.095 
that has an established red sequence. We identified 43 cluster (< z >= 2.095 ±0.004) and 74 field galaxies 
(< z >= 2.195 ±0.083) for which we can measure metallicities. For the first time, we show that there is no 
discernible difference between the mass-metallicity relation of field and cluster galaxies to within 0.02dex. 

Both our field and cluster galaxy mass-metallicity relations are consistent with recent field galaxy studies at 
Z ~ 2. We present hydrodynamical simulations for which we derive mass-metallicity relations for field and 
cluster galaxies. We find at most a O.ldex offset towards more metal-rich simulated cluster galaxies. Our 
results from both simulations and observations are suggestive that environmental effects, if present, are small 
and are secondary to the ongoing inflow and outflow processes that are governed by galaxy halo mass. 

Subject headings: cosmology: observations — galaxies: abundances — galaxies: fundamental parameters — 
galaxies: high-redshift — galaxies: evolution 


1. INTRODUCTION 

The evolution of galaxies is likely dependent on halo mass 
and regulated by baryonic feedback cycles (e.s. JDekel et all 
l2009i ). In addition, environment plays a role in removing and 
truncating gas reservoirs while increasing the frequency of 
galaxy interactions/mergers, which could cause differences in 
morphology and color seen locally for field and cluster galax¬ 
ies (e.g., [Dressier! 1198(1) . However, selecting low redshift 
galaxies by their stellar ages and stellar masses shows that 
their observable properties are independent of environment 
(see iBlanton & Moustakasll2009L for a review). This is also 
confirmed by studies showing that the mass-metallicity rela¬ 
tion of local field and cluster star-forming galaxies are consis- 


tent ( 

Mouhcine et al.,;2007; Ellison et al. 2009; Scudder et al j 

2012 

; Hushes et al. 20131). If there is a difference in the mass- 


metallicity relation as a function of environment at z ~ 0, 
then it is observationally constrained to be around 0.05dex 
(Ellison et al. 2009) while simulations constrain it to be less 


than ~0.05dex (iDave et alj|201 ll) . 

Althou gh it is clear tha t a mass-metallicity relation exists 
at z ~ 2 dErb et ahl 120061 : ISanders et all 12014 ISteidel etakl 
120141k it is unclear if it differs as a function of environment 
at high redshift where clusters are beginning to form. Ex¬ 
ploring the environmental influence on galaxies at high red¬ 
shift is difficult since only a handful of spectroscopically 
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confirme d clus ters with z >, 2 exist (e.g., iKurk et alJ 12004 


Galametz et al. 2013 

Gobat etal.1 2013^ Shimakawa et alJ 

20 3 Valentino et al 

2014; Yuanetal. 2014}). Onlv three 


studies have attempted to address the effects of environment 
on the mass-metallicity relation using small samples of galax¬ 
ies, or inhomgenous field and cluster selections, and have pro¬ 
duced conflicting results. A star-forming protocluster con¬ 
taining 6 galaxies at z = 1.99 was found to be up to 0.25dex 
poore r in metals than their field counterparts (IValentino et all 
120141) . while two star-forming protoclusters at z = 2.2 (24 
galaxies) and at z = 2.5 (36 galaxies) were found to be more 
enriched than those of field galaxies for M* < icy m.t, 

( Shimakawa et all 12014ft . On the other-hand, iKulas et all 
( 20131) found for 20 star-forming protocluster galaxies that 
they do not exhibit a dependence of metallicity on mass (i.e., 
zero slope) with the low-mass protocluster galaxies showing 
an enhancement in metallicity compared to field galaxies. Al¬ 
though there appears to be some disagreement between stud¬ 
ies, these ~0.15dex differences are marginally detected at less 
than two sigma. 

Here, we present Keck/MOSFIRE observations from the 
ZFIRE surve\Q (T. Nanayakkara et al. in prep) in order 
to clarify the debate on the dependence between the mass- 
metallicity relation and environment found at z ~ 2. We iden¬ 
tified 43 cluster (< z >= 2.095 ± 0.0 04) and 74 field ga laxies 
(<z >=2.195 ±0.083) dSnitler et al]l20TllYiian et aMoTi 
for which we can measure metallicities. This study is unique 
given the large number of galaxies obtained for a single clus¬ 
ter, with field galaxies selected the same way and observed si¬ 
multaneously with the same instrument. In §|2]we describe the 
data and our metallicity measurements. In §0 we present the 
mass-metallicity relation for field and cluster galaxies at z ~ 2. 
Our data are consistent with previously published field galaxy 
relations, however contrary to other cluster galaxy works, we 
do not find a difference in the relation as a function of envi¬ 
ronment. We further show that this result is consistent with 

8 http://zfire.tamu.edu 
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Fig. 1.— Representative sample of flux-calibrated ID and 2D MOSFIRE spectra of four cluster members in order of decreasing stellar mass. The spectra are 
plotted with respect to the systemic velocity of each galaxy, as definded by their nebular emission-lines. Both the [NII] doublet and Ho are shown in the ID and 
2D spectra. Strong sky-lines are seen as long virtical lines in the 2D spectra. The ID data (blue). 1 cr error spectra (green) and best-fit Gaussian models (red) 
are shown. The metallicity for each object is also lis ted. The top panels show 4" X 4" three color HST images, usi ng the F814W, F125W and F160W filters, 
obtained from the publicly available CANDELS data iGrogin et al. 2011; Koekemoer et al. 2011; Skelton et al. 2014). 


current simulation predictions. We end with concluding re¬ 
marks in § 0] 

Throughout, we adopt a h = 0.70, 17m = 0.3, = 0.7 cos¬ 

mology. 

2. MOSFIRE SPECTROSCOPIC OBSERVATIONS AND SAMPLE 

The c = 2.095 cluster was previously identified using the 
photometric redshifts (ISpitler et al.1120121) from the medium- 
band FourStar Galaxy Evolution Survey (ZFOURGE, Straat- 
man et al. in prep). As part of the spectroscopic followup 
to ZFOURGE, the ZFIRE survey used Keck/MOSFIRE to 
spectroscopically confirm the existence of a Virgo-like pro¬ 
genitor at z = 2.095 containing at least 57 members with a 
velocity dispersion of 550 km s -1 (lYuan et al.l 120141) . Fur¬ 
thermore, the ZFIRE survey yielded an additional 123 spec¬ 
troscopic redshifts of field galaxies with 1.98 < z < 3.26. 
The spectroscopic targets were selected using the photometric 
redshifts from the K-band selected catalog from ZFOURGE 
(Snitl er et al.l 1201 41 Straatman et al. in prep), which have 
an accuracy of Az/(1 + zspec) = 2% ClTomczak et al.l l2014t 
lYuan et aI1l2014l) . 

The MOSFIRE near-infrared K- and H-band spectroscopic 
observations, data reduction and flux calibration procedures 
are described in lYuan et akl (12014l) with additional details to 


be presented T. Nanayakkara et al. (in prep). All spectra are 
calibrated to vacuum wavelengths. Our typical 3er flux limit 
is 1.8 x 10 -18 ergs/s/cm 2 . 

Gaussian profiles were simultaneously fit to Ha and [NII] 
emission-lines to determine their total flux. The line centers 
and velocity widths were tied together for a given pair of lines. 
Examples of our 2D and ID MOSFIRE spectra and line-fits 
are shown in Figure Q] for galaxies that have stellar musses 
ranging between log(M/M©)=10.4-9.2 (also see lYuan et al.l 
(120141) for additional examples). 

We compute a gas-phase oxygen _abundance for each 
galaxy using the N2 relation of lPettini & Patiell ( 2004 ) where 
12+log(Q/H)=8 90+0. 57 xN2 (N2=log([Nn]/Ha)). The 
iPettini & Page! (120041) relation is established from a sample 
of low-redshift extragalatic H II regions for which both N2 and 
oxygen abundances are directly measured. Applying this re¬ 
lation to high redshift galaxies implies the caveat/assumption 
that the ionization parameters and electron densities are simi¬ 
lar to those of low redshift galaxies, however, it is known that 
typical ionization parameter s are much higher for z ~ 2 galax- 
ies w ith log(LQ>7.3 (e.g ., IShirazi et al.ll2014t ISteidel et al.l 
1201 4t iKe wlev et all 1201 51) . It is important to note however, 
that lKewlev et al.l (120151) has shown that our z ~ 2 field and 
cluster galaxies have consistent ISM conditions, which al- 
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Fig. 2.— (left) The mass-metallicity relation for star-forming z = 2.1 cluster (red) and z ~ 2.1 field galaxies (black). Circles indicate metallicity measurements 
where both Ha and [Nll]are detected with greater than 3er significance and downward arrows are quoted as 1 <t limits on [Nil], The average error in deter¬ 
mining the mass of galaxies using ZFOURGE photometry and FAST of ±0.076 is shown in the botto m left. The green circles show the individual MOSDEF 
detections I Sanders et al. 2014), while the grey shaded region shows the fit to KBSS-MOSFIRE (Steidel et al. 2014). The dotted horizontal line solar abundance 
i Asnlund et alj|200^ i. — (right) The dotted lines are bootstrap fits, with la limits, to the mass-metallicity data from the left panel for cluster (red) and field 
(black) galaxies (see text for details). The solid points are stacked spectra within a given mass bin and note that these are consistent with the fitted data. We 
further show additional binned data from the literature. Note both cluster and field galaxies are consistent within the 1 a errors. 


lows for a meaningful direct comparison between the mass- 
metallicity relations of field and cluster galaxies. 

From the parent sample of lYuan et al.l (120141) . we identify 
49 cluster and 86 field galaxies with an Her flux detected at 
greater than 3er significance. Using the ZFOURGE AGN cat¬ 
alog (M. Cowley et al. in prep), we remove AGNs from our 
sample (7 from the field, 3 from the cluster). An AGN is 
identified as, or acombination of, an infrared source (follow¬ 
ing methods of lDonlev et 41.1^0721), an X-ray source (follow¬ 
ing methods of ISzokolv et al.ll2004l) or a radio source (fol¬ 
lowing methods of G. Rees et al. 2015, in prep). In addition, 
we further require log([N Il]/Ha)< -0.3 JSanders et al.l2014b . 
which removes 8 additional AGN (5 field, 3 cluster). 

Our final sample contains 43 cluster (< z >= 2.095 ±0.004) 
and 74 field galaxies (1.98< z <2.56, < z >= 2.195 ±0.083) 
for which we can measure metallicities. We require a 3cr de¬ 
tection significance level for [Nil], otherwise lcr detection 
limits are shown. Our final sample contains 22 metallicity 
measurements and 21 limits for the cluster galaxies and 41 
metallicity measurements and 33 limits for the field galaxies. 

We use stella r masses computed from the ZFOURGE pho¬ 
tometry using iBruzual & Chario t (120031) stellar population 
models with FAST ( Kriek et all 12009 ). assuming exponen¬ 
tially declining star formation histories, solar metallicity, a 
Chabrier (2003) initial ma ss function and cons trained to the 
spectroscopic redshift (see lTomczak et alll2014l . for details). 

2D Kolmogorov-Smirnov test shows that the cluster and 
field galaxies have statistically consistent (P(KS)< 1 cr) obser¬ 
vational properties such as extinction, (determined by FAST) 
and star-formation rates (determined by Ha - to be discussed 
in Kacprzak et al., in prep). 

3. RESULTS 


3.1. ZFIRE Observations 

In Figure[2](left), we show the mass-metallicity relation for 
our field (black) and cluster (red) galaxies at z ~ 2, including 
la limits. We have metallicity measurements for a signifi¬ 
cant range of galaxy masses from 8.9 <log(M/Mg)< 11.0. 
We further show the 53 individual detection from MOS¬ 
DEF (ISanders et al.ll2014l) and the fitted relation from KBSS- 
MOSFIRE "|Steidd^et^]i20l3). Note that the scatter in the 
data is similar to the scatter in the MOSDEF distribution and 
our data seem to follow the fitted data from KBSS-MOSFIRE. 
Furthermore, our cluster and field galaxies appear to have 
similar mass-metallicity distributions. 

In Figure [2] (right), we show a bootstrap fit (1000 times) 
to the data from Figure [2] (left) using 12+log(0/H)= y, + 
m,(M—10). We fit the relation including the la limits using 
the expectation-maximization maximum-likelihood method 
of [Wolvnetz I (119791) . Including limits we find for the cluster 
y c = 8.370 ± 0.030 and m, = 0.424 ± 0.068 while for the field 
we find yf = 8.384 ± 0.018 and mj = 0.245 ± 0.044. These 
fits, along with their lcr errors, are shown as dotted lines in 
Figure [2] (right). From the fitted data, we find that both clus¬ 
ter and field galaxies exhibit similar mass-metallicity relations 
and are consistent within lcr. The maximum separation be¬ 
tween the fitted field and cluster galaxy metallicity is at most 
0.014±0.035dex. There could be a difference in the slope be¬ 
tween field and cluster galaxies, which is mostly driven by the 
larger number of low mass cluster metallicity limits. 

To explore the similarities between the field and cluster 
galaxy data, we further stack the spectra in five mass bins 
with roughly equal numbers of galaxies per bin. We stacked 
the spectra, weighting by the uncertainty spectrum, to deter¬ 
mine the typical metallicity in a given mass bin shown in Fig- 
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ure |2 We found no discernible difference if we weighted 
each spectrum by its Ha luminosity. Again, the stacked spec¬ 
tra show that cluster and field galaxies have equivalent mass- 
metallicity relations and are consistent within the ler errors. 
The stacked spectra are also consistent with the fits to the in- 
dvidual datapoints. 

The fits and stacked data for both our field and cluster 
galaxi es are consistent w ith the previous resul ts of MOS- 
DEF (ISanders et all Imi KBSS-MOSFIRE (ISteidel et al.1 
120 1 41) and lErb et ak (120061) . However, the low mass end of 
the field galaxies marginally deviate away from the KBSS- 
MOSFIRE fit, but is consistent with their binned data (not 
shown). Any observed minor differences could due to sam¬ 
ple selection biases (KBSS-MOSFIRE is rest-frame ultravio¬ 
let selected while MOSDEF and ZFIRE are rest-frame opti¬ 
cal selected) or how detection significance levels and fits were 
conducted. However, given that our cluster and field data are 
consistent with all previous works, further validates the lack 
of a difference between field and cluster galaxies. 

3.2. Cosmological Simulations 

To further examine if the environment of galaxies actually 
influences their metallicities, we perform a similar analysis 
using cosmological simulations. The simulations are per¬ 
formed with a Gadget-3 based hydrodynamic al code that in¬ 
cludes important baryon physics such as star formation, feed¬ 
back from supernovae and AGN, and chemical enrichment 
from Type II and la supernovae, and asym ptotic giant branch 
stars with Kroupa IMF (see iKob avashi. Snrin gel & Whitel 
[20071 : Taylor & Kobavashl 12014 for the details). Metallic¬ 
ities derived for Kroupa IMF are virtually identical for to 
those derived for a Chabrier IMF. The input nucleosynthesis 
yields are in excellent agreement with the observed elemen¬ 
tal abundances in the Milky Way Galaxy from carbon to zinc 
(lKobavashi& Nakasatoll201 ll). 

Different from the IDave et all (1201 ll) simulations, we in¬ 
clude AGN feedback, whereby AGN-driven winds eject met¬ 
als from massive galaxies to the intergalactic medium, and 
some less-massive galaxies undergo external enrichment de¬ 
pending on the environment. However we find that for the 
mass-metallicity relation, the effect of AGN is small and 
the trend originate mainly from supernova/hypernova-driven 
gal actic winds, which is much larger in less-massive galaxies 
(IKobavashi. Springel & White ! 2007 ). 

The simulations are run for a 25/j -1 Mpc box with two 
different initial conditions. The cluster simulation is chosen 
from ten realization to have a strongest central concentration, 
which gives the most massive galaxy with the stellar mass 
~ 1O 12 M 0 dTavlor & Kobavashil2014l) at z = 0, while the field 
simulation is same as used in IKobavashi. Springel & Whitel 
( 2007 ) and does not contain any central concentrations. The 
stellar masses are estimated by fitting a core-Sersic pro¬ 
file, and the oxygen abundances are measured in 15/r 1 kpc 
weighted by the star formation rates of gas particles to be 
comparable to our emission-line observations. At z = 0, the 
mass-metallicity relations of simulated galaxies are in good 
agreement with the observed mass-metallicity relations both 
for stellar and gas-phase metallicities dTavlor & Kobavashil 
[ 20151) . These relations evolve as a function of time, with 
lower metallicities and a steeper slope at higher redshifts 
(Taylor & Kobayashi 2015, in prep.). 

Figure [3] shows the mass-metallicity relations of the 
“cluster” and “field” simulations at z = 2. We apply an 



Fig. 3.— Mass-metallicity relation for simulated field (black) and cluster 
(red) galaxies. The solid and dashed lines show the moving median with 68% 
of confidence level. The dotted line is an observational limit of our ZFIRE 
survey (see text for details). We find that both field and cluster galaxies have 
consistent mass-metallicity relations with maximum separations of 0.1 dex. 

observationally-based detection limit on the simulated mass- 
metallicity distribution. The N2-based metallicity measure¬ 
ment is confined by the detection limit of [Nil], Based 
on a few simple assumptions, [Yuan et akl ( 2013 ) used the 
SFR-mass relation, flux detection limit and redshift to de¬ 
termine observational detection limit on the mass-metallicity 
relation. For the ZFIRE survey, our typical 3a flux limit 
is 1.8x l(r ls ergs/s/cm 2 resulting in the limiting relation of: 
12+log(0/H) > -O.274xlog(M*/M 0 )+lO.371. This obser¬ 
vational detection limit is important for understanding the 
incompleteness and biases due to observational limitations. 
While all of our observational data reside above this limit, we 
apply it to our simulated data (dotted line). 

In Figure[3] the solid and dashed lines show the moving me¬ 
dian with 68% confidence levels also shown. Both the field 
and cluster galaxies exhibit similar mass-metallicity relation¬ 
ships within la. The maximum separation between the fitted 
field and cluster galaxy metallicity is at most 0.1 dex. This 
is consistent with previous results from z = 0 simulations that 
constrain the difference to be less than ~0.05dex (IDave et all 

HElD)- 

Note that there is an ~ 0.5dex offset in the zero-point be¬ 
tween simulations and observations, which can be partly at¬ 
tributed to t he uncertain normalization in metallicity m easure- 
ments (e.g jKewlev & Ellisonll2008l: IDave et alJl20lH) . 

4. CONCLUSIONS 

High redshift protoclusters are ideal laboratories to study 
the possible influences of environment on galaxies. It is im¬ 
portant to establish whether the environment can affect chem¬ 
ical evolution via restricting inflow and outflow, which may 
be responsible for galaxies migrating from the blue to the red 
sequence, or if a galaxys’ chemical evolution is established 
early and is solely dependent on its internal evolutionary pro- 
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cesses. 

The environment may already be affecting the morpholog¬ 
ical properties of galaxies at z = 2. In our cluster, quies¬ 
cent galaxies have similar colors and sizes relative to field 
quiescent galaxies, however, cluster star-forming galaxies 
are larger and redder on average when compared to field 
star-forming galaxies - suggesting environment can trans¬ 
form j>ome galaxy properties at this early stage of formation 
(1 Allen et al.il20151) . 

It is a concern when searching for metallicity differences as 
a function of environment that the precise metallicities de¬ 
rived for high redshift galaxies are likely incorrect due to 
their different ISM conditions relative to z = 0 galaxies - 
where emission-line metallicity indicators are calibrated (e.g., 
iPettini & Pagell 120041) . However, this does not affect the 
relative offsets between the field and cluster galaxies given 
iKewlev et al. ( 2015 ) has shown that for our cluster, both field 
and cluster galaxies have similar ISM conditions. 

Here, we present the chemical abundance of a Virgo-like 
progenitor at z = 2.1 (lYuan et al.l 120141) that has an estab¬ 
lished red sequence (iSpitler et al.l 120121) . Using our sample 
of field galaxies at a similar redshift, we show that the mass- 
metallicity relation of field and cluster galaxies are consistent 
with lcr errors. We further show that both our field and clus¬ 
ter galaxies have consistent mass-metallicity relations when 
compared to other field galaxy surveys at z ~ 2 (lErb et alJ 
l2006t [Sanders et aPl2014t ISteidel et al.ll2014l) . Although we 
cannot rule out weak environmental trends in chemical en¬ 
richment, our analysis shows that the difference between field 
and cluster galaxies in the mass-metallicity relation is less 
than 0.02dex. 

Our simulations have shown that cluster galaxies may be 
marginally more metal rich than field galaxies by at most 
O.ldex (although consistent within the scatter). This could be 
due to the removal of gas from the outskirts of galaxies which 
can produce an observed increase in metallicity by ^O.ldex 
dHughes et al .1120 1 3i) or due to the shorter gas recycling times 
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Our results from the simulations and observations are sug¬ 
gestive that environmental effects, if present, are secondary to 
the ongoing internal processes within z ~ 2 galaxies that are 
likely governed by halo mass. 
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